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Abstract. Only a small fraction of galactic nuclei in the 
nearby universe are luminous; most of them are currently 
dim. We investigate the feasibility of constraining the X- 
ray luminosity in the recent past (up to ~ 10 6 years ago) 
of the nucleus of a cluster dominant galaxy by measuring 
the contribution of scattered radiation from the central 
source to the surface brightness of the intracluster gas 
dominated by thermal emission. We show that resonance 
X-ray lines present an advantage over the adjacent con- 
tinuum, because the relative contribution of the scattered 
component is typically larger in the line case by a sig- 
nificant factor of 3-10. As an example, we estimate the 
level of constraint that could be derived from future fine 
spectroscopic observations on the past X-ray luminosity 
of the nearby M87 and Cygnus A active galaxies. For in- 
stance, scattered line radiation should be detectable from 
the Virgo cluster if the X-ray luminosity of M87 was a 
few times 10 44 erg s" 1 , or 10~ 3 of its Eddington lumi- 
nosity, until a few times 10 5 years ago. For comparison, 
upper limits inferred from the available XMM-Newton and 
Chandra X-ray continuum data are typically 1-2 orders of 
magnitude weaker. 

The same method can be applied to distant power- 
ful quasars (at redshifts z 1) if they have cluster- 
like gaseous coronae, as suggested by Rosat and Chan- 
dra observations of active galaxies at z <J 1. Their sur- 
face brightness profiles in the X-ray continuum above 
~ 10 kcV ^> kT/(l + z) (where T is the gas tempera- 
ture) should be dominated by rcdshiftcd scattered radia- 
tion from the quasar. Therefore, measurements with forth- 
coming mirror-optics hard X-ray telescopes could give in- 
formation on the lifetime of quasars and parameters of the 
hot gas around them. 

Key words: galaxies: active - quasars: general - scatter- 
ing - X-rays: galaxies: clusters 



1. Introduction 



It is commonly believed that active galactic nuclei (AGN) 
are powered by accretion onto a supermassive black hole. 
However, a detailed description of this process is still miss- 
ing. One of the principal obstacles to making progress here 
is that we know very little from observations about the 
variability of AGN on timescales longer than a few tens of 
years. 

Nuclear variability on the longest timescales up to the 
Hubble time is expected to be governed by the supply of 
gas from galactic scales to the central kpc. Such transport 
can be triggered by mergers and interactions of galaxies 
(e.g. Hernquist 1989| ). There are reasons to believe that 
AGN can be strongly variable on timescales much shorter 
than the Hubble time. For example, Ciotti & Ostriker 

Tabor (|99§), 
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( |1997| , |200lD , in a follow-up to Binney 
have suggested a mechanism for elliptical galaxies in which 
an accreting central black hole heats the ambient gas to 
the point at which the accretion stops, then fueling re- 
sumes after the gas has cooled. In this model, strong out- 
bursts during which the AGN luminosity approaches the 
Eddington critical value can occur every 10 8 — 10 9 years. 
A similar outburst behaviour but on shorter timescales of 
typically 10 6 years can result from thermal-viscous insta- 
bility e xpect ed to operate in AGN accre tion disks (Lin & 
Shield |l98| Mineshi ge fc Shields |l990| : Siemiginiowska, 
Czcrn y & Kostyunin 199€ ; Burderi, King & Szuszkiewicz 
1998 ). We may also mention tidal disruption of stars by 
supermassive black holes, which can give rise to rare (ev- 
ery 10 4 -10 5 years per galaxy) , short (mont hs to years) and 
bright (sub-Eddington) flares (Rees 1988 ). 

Since we are unable to explore directly the long-term 
variability of galactic nuclei, some indirect methods are in 
need. A feasible way would be to observe radiation which 
was emitted in the past by a galactic nucleus and later 
scattered toward us by the interstellar medium. Since the 
characteristic size of a galaxy is a few tens kpc, we could 
see the "echo" of previous AGN activity ~ 10 5 years af- 
ter the central source turned off. This idea has been put 
forward by Sunyaev, Markevitch & Pavlinsky ( 1993 ) in 
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relation to our Galaxy. It was suggested that the diffuse, 
hard X-ray emission observed today in the direction of gi- 
ant molecular complexes in the central 100 pc of the Milky 
Way is scattered radiation emitted by Sgr A* in the past. 
This hypothesis has been confirmed by the detection with 
the ASCA and Chandra satellites of diffuse emission in 
the 6.4 keV fluoresence iron line from the Sgr B2 cloud 
(Koyama et al. 1996; Murakami et al. 2001). Cramphorn 



& Sunyaev ( 2002[ ) have recently elaborated on this prob 



lem and derived upper limits on the X-ray luminosity of 
Sgr A* during the past 10 5 years. 

It seems natural to extend the above approach to ellip- 
tical galaxies and in particular to cluster dominant galax- 
ies. Recent observations with the Hubble Space Telescope 
and ground-based telescopes have confirmed previous sug- 
gestions that 1) powerful low-redshift (z < 0.5) AGN are 
found in luminous and preferentially early-type galaxies 



(Bah call et al. |1997| ; Hoo per, I mpey & Foltz [1997) ; Boyce 
et al |1997[ McLure et al. |1999| ; Schade, Boyle & Letawsky 
2000] ) and 2) they are generally located in cluster environ- 
ments (McLure & Dunlop |2001| ) . Since cluster cores typi- 
cally have sizes r c ~ 100 kpc, we could hope to detect ra- 
diation emitted by the nucleus (plus possibly by a jet) of a 
cluster dominant galaxy until a few times r c /c ~ 10 6 years 
ago and scattered by the hot intracluster gas toward us. 
Up to ~ 1% of the radiation of the central source can be 
Thomson scattered by free electrons within the gas. 

Previous theoretical efforts in this direction have fo- 
cused on the prospects for detecting scattered (strongly 
polarized) AGN radiation at radio and optical wave- 

Gilfanov, Sunyaev 
198C| : 



sions. Furthermore, given the huge collecting area (several 
square meters) of these planned telescopes in the stan- 
dard X-ray band, it should become possible to search 
for scattered X-rays from distant quasars (at redshifts 
z > 1). Indeed, since the radiation from a quasar gets 
redshifted on its way to us, a significant gain in the ra- 
tio F AGN /F gas cx [E(l + z)]- a exp[£(l + z)/kT] could be 
achieved already at E ~ 10 keV. Such observations would 
permit one to constrain the lifetime of quasars and the 
parameters of the hot gas around them. This possibility 
is one of the issues discussed in the present paper (partic- 
ularly in §P). 

On the observational side, evidence has been accumu- 
lated over recent years that AGN are often surrounded 
by gas atmospheres typical of clusters of galaxies. Early 
ROSAT observations indicated the presence of luminous 
extended X-ray emission around several powerful quasars 
with redshifts out to z = 0.73 (Crawford et al. 1999 ; Hard- 
castle & Worrall 1999). More recently, the Chandra ob- 
servatory clearly resolved extended hot gas around sev- 
eral relatively nearby AGN, including 3C295 (z = 0.461, 
Har ris et al. |2000| ), 3C220.1 (z = 0.62, Worral et al. 
200l |200lD, 3C 123 (z = 0.2177, Hardcastle, Bir kinsh aw & 
Worrall |200iT ), Cyg A (z = 0.0562, Smith et al. [2002D , and 



2002; ). 



H1821+643 (z = 0.297, Fang et al. 

Chandra has also detected luminous extended (to radii 
~ 100 kpc) X-ray emission surrounding two high-redshift 
powerful radio galaxies 3C294 at z = 1.786 (Fabian et 



lengths (Sunyaev 1982 



1987a; Sholomitskii & Yaskovich 



Churazov 
Wise & Sarazin 



al. |2001| ) and PKS 1138-262 at z = 2.156 (Carilly et al. 
|2002 ). Assuming for these objects that the emission of the 
central source is isotropic and its X-ray luminosity has 
been the same in the past as it is now (~ a few 10 45 erg 



1990|7|1992|; Sarazin & Wise Il992; Murphy & Chernoff s x ), we are able to estimate (see gT2[) that the gas surface 



1993| ). However, it would be particularly interesting to get 



information on the history of the central source from ob- 
servations of X-ray scattered radiation, since the bulk of 
the AGN X-ray emission is believed to originate within a 
few gravitational radii of the central black hole. 

One feasible way to detect the scattered X-rays against 
the thermal emission of the hot gas would be to perform 
observations at £ > kT, where E is the photon energy 
and kT ~ 1-10 keV is the gas temperature. Indeed, the 
AGN spectrum is a power-law, Fagn oc E~ a with a ~ 1, 
while the plasma bremsstrahlung spectrum is exponen- 
tially declining, F gas cx exp(-_E/fcT), at E > kT. Unfor- 
tunately, the angular resolution and sensitivity of present- 
day hard X-ray telescopes are not sufficient for mapping 
clusters of galaxies and detecting scattered AGN radia- 
tion. However, the situation may improve dramatically 
with the advent of grazing incidence X-ray telescopes sen- 
sitive up to ~ 40 keV, such as those under considera- 
tion for the projected Constellation- x|] and XEUSQ mis- 



brightness should be dominated by scattered power-law 
radiation from the AGN already at E ;> 10 keV. 

In this paper we also discuss another possibility of de- 
tecting scattered AGN X-ray radiation. We show (in ^ 
and §[3^) that at least for nearby AGN located in cluster 
dominant galaxies or in isolated giant elliptical galaxies, 
one can gain significantly by measuring the scattered sur- 
face brightness in resonance X-ray lines rather than in the 
continuum. Such observations would ideally require high 
(~ a few eV) spectral resolution in order to separate the 
resonance lines from other lines. The resolution needed to- 
gether with high sensitivity can be achieved with future 
high-energy astrophysics observatories such as Astro-E2Q, 
Constcllation-X and XEUS. As an example, we estimate 
(in §^) the level of constraints it should be possible to 
derive on the past X-ray luminosity of the nearby M87 
and Cyg A active galaxies. Yet better constraints could 
be obtained by measuring the polarization of the scat- 
tered X-ray resonance line radiation (see Sazonov et 
al. |2002D . 



1 http:/ /constellation. gsfc.nasa.gov/docs/main. html 

2 http:/ /astro. esa.int/SA-general/Projects/XEUS/main/ 
xeus_main.html 



http:/ /www.astro. isas.ac.jp/astroe/index-e.html 
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2. Advantage of scattered X-ray lines over the 
scattered continuum 

Consider an AGN surrounded by hot (T ^ 10 7 K) gas. 
The following emission components will contribute to the 
X-ray surface brightness of the gas: 

— First, there is X-ray emission (line plus continuum) 
of the hot plasma in which collisional ionization and 
excitation dominate. However, a bright AGN may ad- 
ditionally photoionize the plasma. In our calculations 
below we fully neglect this effect, although it can 
be important when the AGN luminosity is very high 
(quasar- like) , as will be discussed in a separate paper 
(Sazonov & Sunyaev, in preparation). 

— Second, there is radiation emitted by the AGN and 
Thomson scattered by free electrons in the gas. This 
scattered radiation contributes to the surface bright- 
ness in the spectral continuum. 

— Finally, AGN photons with energies falling within the 
cores of X-ray lines of ions of heavy elements can res- 
onantly scatter on these ions, thus contributing to the 
surface brightness in the lines. 

Consider an optically thin volume of gas with elec- 
tron temperature T and number density n , exposed to 
an external spectral flux Fe (E) (measured in units of erg 
cm~ 2 s _1 keV -1 ). 

The plasma continuum spectral emissivity due to 
bremsstrahlung is given by (e.g. Zombeck 1990 ) 

e B ,co„t(£) = 2.3 x IQ-^T-^expi-E/kTjnlgBiT^) 



(erg cm 



''s^keV- 1 ),^) 



where g^ is the corresponding Gaunt factor. 

The energy-integrated emissivity in a resonance line 
by ions of type z due to electron collisional excitation of 
an electron from the ground level i to an excited level k is 
given by (e.g. Zombeck |1990| ) 



en™ = 2.7 x 10- 15 T- 1 / 2 exp(~^ ifc /fcT) 
x n c n z (T)fag ik (T) (ergcm~ 3 s" 1 ), 



(2) 



where n z is the ion number density, which depends on 
the gas temperature. The remaining quantities appearing 
in equation (^) , namely E^. , fa and go- , characterize the 
line itself, denoting its energy, oscillator strength and ex- 
citation Gaunt factor, respectively. 

The rate (per unit volume) of Thomson scattering of 
the external radiation, integrated over the scattering angle 
0, is 



-.E.cont 



(E) 



| j d cos 0(1 + cos 2 9)<%$ >caat (E, I 
= F E (E)n c a T (erg cm -3 s" 1 keV -1 ) 



(3) 



where or = 6.65 x 10 25 is the Thomson scattering cross 
section. Note that we ignore throughout the Klein-Nishina 



correction to the scattering cross section as well as the 
change in the photon energy by scattering, which is justi- 
fiable if E -C mc 2 and kT <C mc 2 , where mc 2 = 511 keV 
is the electron rest energy. 

The energy-integrated rate of resonant scattering (in- 
tegrated over the scattering angle) of the external radia- 
tion in the line i — > k by ions z is 



4.1 x 10" i8 7r 



x F E (E lk )n z (T)r c cfa (erg cm 3 s x ), 



(4) 



where r c — 2.82 x 10~ 13 cm is the classical electron ra- 
dius, c is the speed of light, and the numerical coefficient 
accounts for the unit transition from Hz to keV. 

From equations ([!]) and (^) we can find the ratio of the 
Thomson scattering rate to the bremsstrahlung emissivity: 



c £,cont 
£S,cont 



(E) = 2.9 x 10 



-5 



x Fe^T 1 ' 2 exv{E/kT)n- l g^{T,E). (5) 

The ratio of the resonant scattering rate to the emis- 
sivity in the line z — s- fc is given by a similar expression, 
which follows from equations (0) and (H), 



line 



4.1 x 10" 



E. 



x F E (E ik )T x l 2 eMEik/kT)n- x g^(T). (6) 
We can now compare the ratios @ and (||) for E = 



cat ^scat 



R = 



— c lino 



£ _E,cont _ -j A 9 B(T,Ejk) 



= 1.4 



Qinc ^i^cont 



(7) 



In the low-temperature limit (kT <gC Ei k ), the excita- 
tion Gaunt factor g ik ranges between 0.1 and 0.25 for most 
interesting lines; it slowly increases with temperature and 
reaches 0.2-0.3 at kT ~ E ik (e.g. Mewe, Gronenschild 



& van der Oord 1985 ). The bremsstrahlung Gaunt factor 
also is a weak function of temperature: gs = (kT/E) 0A 
when E ~ kT. Summarizing these facts, we may simply 
write 



R - 6 for kT - E, 



ik • 



(8) 



This approximate relation is accurate to within a factor 
of 3 for all resonance X-ray lines of interest to us. Fig. [I] 
shows, for several resonance lines, the dependence R(T) 
computed using the MEKAL plasma emission code in 
combination with other atomic data sources (see §|] for 
details). In these computations the contribution to eii no 
of unresolved (defined as those with energies within one 
Doppler width of the resonance energy) satellite lines was 
taken into account. This effectively leads to a smaller value 
of R compared to the definition (fjj). However, this reduc- 
tion is only significant for the 6.70 keV and 6.97 keV Ka 
lines or iron when kT -c Ei k . We should note that the 
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Fig. 1. Ratio R defined by eq. [[7| as a function of plasma 
temperature for various resonance X-ray lines: Fe XXIII 
L-shell line at 1.13 keV (solid), Fe XXIV L-shell line at 
1.17 keV (dotted), Si XIII Ka line at 1.86 keV (short- 
dashed), S XV Ka line at 2.46 keV (long-dashed), Fe XXV 
Ka line at 6.70 keV (dot-short-dashed), Fe XXVI Lya line 
at 6.97 keV (dot-long-dashed), and Fe XXV K/3 line at 
7.88 keV (short-dash-long-dashed). 

large values (R ~ 20 at kT ~ Eik) obtained for the iron 
L-shell lines at 1.13 keV and 1.17 keV may be partly the 
result of the different values for the oscillation strengths 
being used in the MEKAL code with which we compute 
eii n e and in the list of resonance lines of Verner, Verner & 



Fcrland (1996) which we use to compute ef™ c . We do not 
attempt in this paper to correct for this and possibly other 
inconsistencies caused by the simultaneous use of several 
sources of atomic data. However, we estimate that some 
of our computational results obtained below (in §|J) may 
contain a relative systematic error of < 2. 

We point out that R is typically several times smaller 
for intercombination lines than for resonance lines and 
R « for forbidden lines. For this reason we consider 
only resonance lines throughout. 

We have thus found that for a given resonance line with 
energy Eik > the relative contribution of scattered external 
radiation to the volume X-ray emissivity of hot gas is typi- 
cally larger by a significant factor of the order of 3-10 than 
for the continuum at Eik- We also know that for typical 
cluster temperatures kT ~ a few keV the combined surface 
brightness of the intracluster gas in the resonance X-ray 
lines with Eik ~ kT is comparable to that in the contin- 



Fig. 2. Sketch of the model. The points with equal travel 
time from the AGN to the observer lie on a paraboloid. In 
the switch-off case, the thermal plasma located in the grey 
and black regions is scattering the radiation emitted by 
the AGN. In the case of a flare of the AGN, only the grey 
area is filled by photons and contributing to the scattered 
emission. A possible path for a photon scattered into the 
line of sight of the observer is marked by a dashed line. 
The system is rotationally symmetric about the z-axis. 

uum. We therefore arrive at the conclusion that resonance 
X-ray lines present significant advantage over the contin- 
uum for searches of scattered AGN radiation in clusters 
of galaxies and elliptical galaxies performed at E ~ kT. 

After we have introduced the R factor in equation ^ 
and shown that this factor is only weakly dependent on the 
resonance line and on gas temperature, we can continue 
our treatment in parallel for continuum and line emission. 

3. AGN in the center of a beta-cluster 

We shall now consider the following model (see Fig. ||) . A 
supermassive black hole is located in the center of a clus- 
ter of galaxies containing hot gas with a b eta-la w radial 
density profile (Cavaliere & Fusco-Femiano 1976 ): 

n c = n (l + r 2 /r c 2 )- 3 ' 3 / 2 , (9) 

where no is the central electron number density and r c 
is the core-radius. The gas temperature T and element 
abundances are constant over the cluster. 



At a moment t a 



-A (as measured by an earth 



bound observer) the central source switches on (i.e. be- 
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comes an AGN) and at a later moment t g = it switches 
off. The AGN emits isotropically X-ray radiation having 
a power-law spectrum of photon index 7, at a constant 
luminosity Lx in the energy band [E\, E^]: 



rE 2 C l2_ 7 

A I E-^EdE = A-± 



E 2 ^ 



7 



(10) 



The spectral flux at a given distance r from the AGN is 
therefore 



F K = 



47r(K 2 - 7 



(2 - j)L x 



E 2 - 7 )r 2 ~ 



(11) 



Our goal is to find the surface brightness profile of 
the cluster, including the contribution of scattered AGN 
radiation, at a given moment of time t > t on . 

3.1. Approximation of isotropic scattering 

We shall use throughout this paper the approximation of 
isotropic scattering. This simplifies the treatment but re- 
quires explanation. We know that in reality Thomson scat- 
tering occurs according to the Rayleigh phase function 
p(cos9) = 3(1 + cos 2 0)/8. In the case of resonant scat- 
tering, the phase function depends on the line and gener- 
ally can be represented as a weighted sum of the Rayleigh 
phase function and the isotro pic function p(cos6>) = 1/2 
(see e.g. Chandrasekhar 1950 ). However, for our given ge- 
ometry of the problem the resulting brightness profiles de- 
pend only weakly (typically by less than 10%) on the phase 
function. For this reason, the approximation of isotropic 
scattering is quite reasonable. 

3.2. Continuum radiation 

We shall first consider the X-ray spectral continuum, for 
which the relevant scattering mechanism is Thomson scat- 
tering. Since the Thomson optical depths of clusters of 
galaxies and elliptical galaxies are small, tt <J 0.01, we can 
use the single-scattering approximation and make some 
analytic estimates. 

The continuum (bremsstrahlung) surface brightness of 
a beta-cluster is described, as a function of projected dis- 
tance p from the nucleus, by the well-known formula, 
which can be obtained by integrating the plasma emis- 
sion along a given line of sight (i.e. along the coordinate 
z - see Fig. ||) using equations ([!]) and (||): 

eE,cont(0 



B 



£\cont 



2.9 x 10 



kT 



■ dz 
Air 

_ 6 r(3/?-i/2) 



r(3/3) 



( n ° V 

V 0.01 cm" 3 J 



100 kpc 



-1/2 



IkeV 

x (l + pi/rl)- 3 ^ 2 

(erg cm~ 2 s _1 sr _1 keV -1 ) 



(E/kT)- { > A ex-p{-E/kT) 



where r = (p 2 + z 2 ) 1 / 2 and we have adopted gs(T,E) = 
(E/kT)-° A . 

The surface brightness of the Thomson scattered AGN 
radiation is given by 



nscat 
-°.E.cont 



^if^ = £/ ft(r)n.(r)«Ml3) 



and our assumption of 



where we have used equation 
isotropic scattering. 

The time-dependent integration limits Z\{t, p) and 
z% (t, p) in equation ( |l3] ) are determined by the loci of scat- 
tering sites giving a fixed time delay t. If the distance be- 
tween the emitter and the observer is much larger than the 
characteristic size of the scattering cloud, which is true in 
our case, the constant-delay surface is a paraboloid with 
its focus at the p osition of the source (see e.g. Sunyaev 
& Churazov 1998 and references therein to earlier work 



considering similar light-echo problems): 



Therefore, 



El 

ct 



ct 



c(t + A) 



c(t + A) 



(p 2 /ct 

OO 



ct)/2 



ifO < t 

if - A < t < 0. 



(14) 



(15) 



(16) 



The integration limits in equation ( |13| ) depend on the 
outburst duration A. We shall first consider two limiting 
cases for which analytic treatment is possible. 

3.2.1. Stationary case 

Let the AGN be a persistent source. The limits of the in- 
tegral e^ont dz are then Zi(t, p) = —00 and 22 (t, p) = 
00, and it is possible, using equation (^|) for n c (r) and 
equation ([ll]) for Fe, to express the integral through hy- 
pergeometric functions. In the particular interesting case 
of j3 = 2/3, the scattered surface brightness is 



Dscat / a 
°E.cont\P 



2/3) 

100 kpc 2-7 
rZ lO 2 -'' - 1 



6.0 x 10' 



L 



x 



M 



BH 



"0 



LwAd 10 9 M m 0.01 cm" 3 



r c /p-(l + p 2 /r 2 )-^ 2 



J5I-7 



(erg cm" 2 s" 1 sr" 1 keV _1 ),(17) 



where Mbh is the mass of the central black hole, L^dd = 
1.4 x 1O 38 Mbh/M is the corresponding Eddington lu- 
minosity, Lx is the AGN luminosity in the energy range 
1-10 keV, and E is measured in keV. 

The ratio of the scattered surface brightness to the 
intrinsic surface brightness of the gas in the continuum is 



oscat 
■"B.cont 



(12) Be" 



COIlt 



(/3 = 2/3) = 2.3 



L x Mbh 0.01cm" 
L^d~d~ 10 9 M Q ~ 
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X(l+/A-c) 



100 kpc 

2 /„2\3/2 



kT \ 01 2-7 
1 keV J lO 2 " 7 - 1 

r c /p-(l + p> c 2 r 1/2 

x£; 1 - 4 - 7 exp(S/fcT) 



(18) 



We can also write down asymptotic expressions appli- 
cable in the limit of p 3> r c for arbitrary values of (3: 



B 



£,cont(P>^c) = 3.2x 10 



_ 5 r(3/3/2 + l/2) L 



x 



AL 



BH 



?? 



T(3/3/2 + l) L Ed dl0 9 M Q 
100 kpc 2-7 ^^^j^j 



oscat 
-£?E con t 



0.01cm" 3 



0> r c ) = 1.1 



10 2 -t - 1 



(p/r c )- 



(ergcm^^-^r^^eV -1 ); (19) 



T(3/3/2 + l/2)r(3/3) 



r(3/3/2 + l)r(3/3 - 1/2) L Edd 

3 2 / j.rri \ 0.1 



M BH 0.01 cm" 3 / 100 kpc 



10 9 M Q 
2 



IkeV 



10 2 ~f - 1 



(p/r c ) 3/3 " 2 S 1 - 4 - 7 exp(£:/fcT). (20) 



It can also be shown that within the cluster core, 
£f,cont (x p- 1 and Bf^t/^cont oc p" 1 for any /?. This 
can be directly verified for (3 — 2/3 using equation ([l8|). 

Fig. ||a demonstrates the dependence of the surface 
brightness profile of a cluster with a central AGN on the (3 
parameter. The scattered surface brightness profiles were 
derived by numerical integration of equation ([13]). Fig. |] 
shows the corresponding plots for the ratio of the scattered 
surface brightness to the intrinsic surface brightness of the 
intracluster gas. We see both from Fig. ^ and equation 
( p(i| ) that outside the core region (at p ;> r c ), the relative 
contribution of scattered AGN radiation to the surface 
brightness increases with p when (3 > 2/3. 

The (3 = 1/3 case presented in Figs. ||a and || cor- 
responds to cluster cooling flows. Typically, gas density 
grows oc r" 1 toward the center of the cooling flow, which 
roughly corresponds to a beta-model with a small r c , large 
n and (3 ~ 1/3. We then find from equations (|l^) and ( |l9| ) 
that i?B jCO nt oc p~ x and S^ont oc p~ 2 (assuming a con- 
stant temperature) for cooling flows. Interestingly, the sur- 
face brightness profile of scattered AGN radiation within 
the core of a beta-cluster (see above) is similar (cx p" 1 ) 
to the intrinsic brightness profile of an (isothermal) cool- 
ing flow. Of course, a cluster core dominated by scattered 
emission from a central powerful quasar will have a very 
different spectrum (power law) than a cooling flow. 

The energy dependent factor E 1A_1 exp(E / kT) ap- 
pearing in equations ( |l8| ) and ( p0| ) has a minimum (for 
7 > 1.4) at E — (7 — lA)kT and increases exponentially 
when E kT. The origin of this dependence is obvi- 
ous: we are comparing the slowly declining (power-law) 
spectrum of the AGN emission with the exponentially de- 
caying efficiency of collisional radiative mechanisms op- 
erating in the intracluster gas. It is clear that for given 
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Fig. 3. a) Surface brightness profiles of scattered AGN ra- 
diation (solid lines) and plasma thermal emission (dotted 
lines) in the X-ray continuum at E = 20 keV, for a model 
in which a central 10 9 M Q black hole is a persistent X-ray 
source with Lx = ^Edd and 7 = 2. The host cluster has 
the following parameters: r c — 100 kpc, n — 0.01 cm" 3 , 
kT = 3 keV, and (from top to bottom) (3 = 1/3, 2/3 
and 1, and is local (z < 1). The dashed lines are p ^S> r c 
asymptotes given by eq. fll9|] , applicable at p > 2r c . In the 
opposite limit of p <C 1 the scattered surface brightness is 
proportional to p" 1 and is independent of f3, as indicated 
in the upper left corner of the figure, b) Surface bright- 
ness profiles of scattered AGN radiation (solid lines) at 
E = 20 keV, measured at different times (t = r c /c, 3r c /c, 
9r c /c and 19r c /c - from top to bottom) after the end of a 
short (A = r c /c) flare with Lx = L^dd- For the same val- 
ues of the parameters as in (a) and (3 = 2/3. The dashed 
lines are p, ct r c asymptotes given by eq. |2^] . The dot- 
ted line is the surface brightness profile of plasma thermal 
emission, c) Same as (b), but for the case where the cen- 
tral source was on (Lx = -^Edd) until t Q ff = (profiles 
order from top to bottom), d) Same as (c), but for the 
case where the central AGN switched on at t = (profiles 
order from bottom to top). 



AGN luminosity and parameters of the gas, a significant 
gain (in terms of the relative contribution of scattered 
AGN radiation to the surface brightness) can be achieved 
by performing observations at E 3> kT, like in our ex- 
amples in Figs. [5] and ^. A cautionary note is necessary 
here. If one wants to constrain the past luminosity of an 
AGN by mapping its host elliptical galaxy in hard X-rays, 
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index for the AGN emission 7 ~ 
depth of the cluster is ~ 0.004. 
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Fig. 4. Radial profiles of the ratio of the surface brightness 
of scattered AGN radiation to the intrinsic brightness of 
the beta-cluster gas, for the same model as in Fig. ||a. The 
solid lines correspond (from bottom to top) to j3 = 1/3, 
2/3 and 1. The dashed lines are the corresponding p^> r c 
asymptotes given by eq. [f20|. 



a restriction will appear due to a significant contribution 
to the measured surface brightness from the (unresolved) 
X-ray binaries of the galaxy. We estimate that in prac- 
tice it should be very difficult to obtain limits better than 
Lx ~ 10 _3 iEdd- Of course, this complication is irrelevant 
for surface brightness measurements carried out over re- 
gions of a galaxy cluster that are almost devoid of point 
X-ray sources, say at distances ~ 100 kpc from the nucleus 
of a central dominant galaxy. 

3.2.2. The case of a high-redshift AGN 

Furthermore, if a powerful quasar is observed at high red- 
shift z, the relevant factor will be [^(l + z)] 1 - 4 ^ 7 exp[E(l+ 
z)/kT]. Consider for example the radio galaxy 3C294 lo- 
cated at z = 1.786, for which a surrounding intracluster 
medium has recently been revealed with Chandra (Fabian 
et al. |200ip . Since the total number of counts registered 
from the source in the 20 ksec Chandra observation was 
small, ~ 100, of which only ~ 30 are from the cen- 
tral AGN, it is only possible to make rough estimates in 
this case. Fabian et al. find that the gas temperature is 
~ 5 keV, the 2-10 keV rest-frame luminosity of the hot 
gas is ~ 2.5 x 10 44 erg s _1 , the corresponding (isotropic) 
luminosity of the AGN is ~ 1.1 x 10 45 erg s^ 1 , the photon 



We can then estimate that the ratio of the flux of scat- 
tered AGN radiation to that of thermal bremsstrahlung 
emission from the gas is only ~ 1% at E <; 2 keV, ~ 3% 
at 5 keV, ~ 30% at 10 keV, - 1 at 12 keV and - 4 at 
15 keV. These estimates assume that the central source 
is isotropic and its luminosity has been the same in the 
recent past (over ~ a few 10 5 years), and are made for the 
total flux from the whole observed X-ray halo (of ~ 15 arc- 
sec, or equivalently ~ 100 kpc radius). We thus see that it 
is desirable to go to energies above 10 keV in searching for 
scattered AGN radiation around the 3C294 galaxy. How- 
ever, photon statistics then becomes a limiting factor. We 
can roughly estimate that the XEUS observatory with its 
effective area of ~ 2 m 2 at 10 keV and angular resolution 
of a few arcsec should be able to detect only a total of ~ 10 
scattered AGN photons above 10 keV in a 100 ks observa- 
tion of 3C294. This, combined with much more accurate 
measurements of the surface brightness of the intraclus- 
ter gas in the standard X-ray band, could be only about 
enough to put constraints on the past history of the 3C294 
central X-ray luminosity. 

Another similar example is the radio galaxy 
PKS 1138—262 at z = 2.2, from which Chandra also has 
observed luminous extended X-ray emission (Carilly et 
al. 2002 ). In this case, the gas X-ray luminosity is sim- 
ilar to that of the 3C294 cluster (unfortunately, the gas 
temperature is unconstrained), but the inferred 2-10 keV 
(rest-frame) luminosity of the central source is four times 
higher, ~ 4 x 10 45 erg s _1 . Moreover, in this case the gas 
emission is likely dominated by scattered AGN radiation 
already at E J> 5 keV (due to the higher redshift and lu- 
minosity), so that a few hundred scattered AGN photons 
could be detected above 5 keV in a 100 ks XEUS-2 obser- 
vation of PKS 1138—262 (taking into account the larger 
collecting area of the telescope at lower energies) . 

We note that photon statistics should be less of a prob- 
lem for more nearby AGN and quasars such as 3C 273 
located at z = 0.158. In fact, it is possible, given the huge 
X-ray luminos ity of 3C 273, ~ 2x 10 46 erg s _1 (e.g. Yaqoob 
& Serlemitsos 1990 ), that a hot intracluster medium sur- 
rounding it, if any, will be first observed in scattered rather 
its own X-ray emission. 



3.2.3. Short AGN flare 

Suppose now that an AGN experienced an outburst in the 
past, between t on — —A and t Q s = 0. The "short" here 
means that the line-of-sight depth of the illuminated (see 
Fig. |2|) zone of the cluster must be much smaller than the 
characteristic size at a given projection radius. Therefore, 
one of the following conditions must be fulfilled: A <C r c /c 
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if p <J r c or A <C p/c if p ;> r c . Then, as follows from 
equations ([l5]) and (|l6|), 

5z(t,p)=z 2 -z 1 ^^(l + p 2 /c 2 t 2 ), (21) 

aild /»? £ E ,cont(P> Z ) rfz » £ f ,cont(P. We tllCU find 

in the limit of p, ct r c 

Bffont = 3.2 X 1(T 5 X 4 3/3/2 L * M BH 

n 100 kpc 2-7 



0.01cm" 



lO 2 ^ - 1 



x (A/t)(r c /ct) 30+1 (l + P 2 /c 2 t 2 )- zfi - 1 E 1 -^ 
(ergcm^s^sr^keV -1 ); (22) 



Dscat 

^g.cont = 12 x 4 3/3/2 



r(3/3) 



BH 



.E.cont 



0.01 cm" 3 / 100 kpc 



r(3/3 - 1/2) L Edd 10 9 M e 
o.i 



kT 



n \ r c J \ 1 keV j 10 2 -t - 1 



7 



x {A/t)(ct/r c ) a P-'(p/ct) 
x £; 1 - 4 - 7 cxp(£;/fcT). 



(23) 



Comparing equations ( J23|) and (|2C|), we see that for 
a given luminosity the contribution of scattered AGN 
radiation to the surface brightness at p ~ ct is smaller by a 
factor of A/i for the flare campared to the stationary case. 
This happens because in the former case the scattered 
radiation comes from a gas layer of depth Sz ~ cA along 
the line of sight, whereas 5z ~ p ~ ct in the latter case. 

Fig. ||b shows a sequence of scattered surface bright- 
ness profiles that would be measured at different times 
after a short AGN outburst. We see that distant regions 
of the cluster reveal no scattered light early on after the 
flare but become progressively brighter as time goes by 
and photons of the flare propagate through the gas. The 
ratio B^^ ont / Be, cont has a maximum at p max ~ ct, as 
indicated by equation (^3|). For example, p max = r c when 
f3 = 2/3. Therefore, the zone of largest contribution of 
scattered AGN radiation propagates outward at an appar- 
ent velocity equal to the speed of light. It is important for 
observations that this zone is broad — 5p c s ~ Pmax ~ 

We also notice that the amplitude of the effect is pro- 
portional to the product L^A, i.e. to the total energy 
emitted by the AGN during the flare. Thus, for example, 
a 10 5 -year outburst at Lx = O.OlLEdd would produce the 
same scattered surface brightness profile (at p, ct ^> r c ) 
as a 10 4 -year flare at Lx = O.liEdd- 

3.2.4. Switch-off and switch-on scenarios 

We next consider a scenario in which the central source 
was persistently (over a time that is longer than the char- 
acteristic light travel time of the cluster) bright in the past 
until it suddenly turned off. A sequence of resulting scat- 
tered surface brightness profiles corresponding to different 



times after the switch-off (t > 0) is shown in Fig. ||c. One 
can see that these profiles differ markedly from those cor- 
responding to the flare case (see Fig. ||b). In particular, as 
time goes by the surface brightness decreases in the central 
region, while it remains practically unchanged in more dis- 
tant parts of the cluster still for a long time ~ p/c after the 
switch-off. This gives rise to a broad (Sp ~ p) maximum of 
the scattered/intrinsic brightness ratio, which propagates 
outward at an apparent speed of light. 

Another interesting situation takes place when the cen- 
tral source is bright at the time of observation but it has 
turned on recently. This "switch-on" case may actually 
correspond to the observed quasars, which may have been 
luminous for a time short compared to the light travel 
time of the gas around them. This situation is depicted in 
Fig. |d. One can see that the radial distribution of scat- 
tered AGN radiation is effectively truncated at a certain 
radius ~ c(t — t on ). This gives the possibility to estimate 
the time during which the quasar has been active. 

3.3. Resonance lines 

The results obtained above can be directly extended to 
the case of resonance lines provided the intracluster gas 
is optically thin (r <§; 1) to resonant sca tterin g. We shall 
consider this limit below and then (in § 3.3.1 ) point out 
differences that may arise when r ^ 1. 

In the absence of external illumination, the equivalent 
width of a line would be 



EW = 



eiinecb/ / e E ,cont(Eik) dz 



1.2 x 10 5 ^/, fc -4^r (keV), 



'g B (T,Ei, 



(24) 



under our assumptions of constant temperature and ele- 
ment abundances. 

Resonant scattering of the AGN photons leads to an 
increase of EW . At the same time, Thomson scattering of 
AGN radiation tends to enhance the continuum near the 
line and thus to decrease its equivalent width. As we know 
from §|^, the first effect is a factor of R ~ 6 larger than the 
second, hence the net effect of the scattering of AGN emis- 
sion in the intracluster gas should be increased equivalent 
widths of the resonance lines. We can then write: 



EW 



CO rz 2 

eiino dz + / ej££ dz ) / 

— oo J Z\ 

e_B,cont dz + I ££,cont dz 
-oo J Z\ 

z 2 ,scat 



EW 



1 + R Q ggont dz/ eg,cont dz 

1 + C £ f a cont dz/ ££,co„t dz ' 



(25) 



If the contribution of scattered AGN radiation to the 
surface brightness of the cluster is small, the fractional 
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increase of the equivalent line width will be 
AEW EW - EW 



EW Q 



EW 

p 1 f z 2 scat 
it — 1 J Zl fc .E,cont 



(E ik ) dz 



5 CE,cont(Eik) dz 



(26) 



Usin g eq uation (pq), one can readily apply all the re- 
sults of § |3.2| to optically thin resonance lines. 

It should be emphasized that equation ( pq 1 ) is strictly 
valid only in the case where both the plasma temperature 
(which affects the ionization balance and plasma emis- 
sivity) and the element abundances are constant along 
the line of sight corresponding to a given projected dis- 
tance p. Indeed, only when these conditions are met, 
eiineW oc e CO nt(r) and ef£f*(r) cx e s c c ^ t {r) ~ see equations 
® _ (q)- ^ n the general case, a coefficient somewhat dif- 
ferent than R — 1 will relate the ratios AEW/ EWq and 

<f,cont dz l I-oo 6 E, cont dz . 

We should also remind the reader that we ignore 
throughout the effect of the AGN on the ionization bal- 
ance (n z ) in the intracluster medium. The irradiation of 
the gas by a very bright central source will first of all alter 
the equivalent width of the emission line relative to the 
value predicted by equation (gjj); it will also affect (how- 
ever to a smaller degree, because both ei; ne and ef-^ are 
proportional to n z ) the ratio AEW/EWq- 

3.3.1. Finite depth effects 

The intracluster gas may be moderately optically thick in 
the stronger X-ray lines (r ^> 1, e.g. Gilfanov, Sunyaev & 
Churazov 1987b), so the analytic results obtained above 
will be only approximately valid in that case. 

The optical depth at the center of a resonance line 
along a given direction I through the cluster is 



tq = / (T Q n z (r) dl- 



VO 



At plasma temperatures typical of galaxy clusters 
(kT ~ 1-10 keV), all interesting X-ray lines have nearly 
Doppler absorption profiles whose width is determined by 
the velocities of thermal and turbulent motions, because 
the line natural width is relatively small. Therefore, the 
cross section at line center is 



CO 



y/nhrecfi. 



AE D 

(h is Planck's constant) 
2kT 



(28) 



with the Doppler width being 



AE D = E lk 



V 2 

"turb 



E, 



AnipC 

2kT 
Am p c' 



1/2 



(1 + 1.4AM 2 ) 



1/2 



(29) 



where A is the element atomic mass, m p is the proton 
mass, Kurb is the characteristic turbulent velocity and M 



is the corresponding Mach number. The role of the ad- 
ditional line broadening due to turbulence becomes more 
important with increasing A. 

Taking into account turbulence line broadening is im- 
portant for two reasons. First, this phenomenon may cause 
some resonance lines, which would otherwise have r > 1, 
to become optically thin. Secondly, some satellite lines 
may fall within the core of the resonance line, thereby 
effectively reducing the factor R defined in equation ((?]). 
As was noted in this effect is particularly important 
for the Ka lines of H-like and He- like iron at kT ~ 1 keV. 

If To is non-negligible, first of all the intrinsic surface 
brightness profile of the intracluster gas will be distorted 
due to diffusion of photons from the central region to the 
outer parts of the cluster. The equivalent line width will 
be somewhat smaller than predicted by equation ( p4| ) at 
p < r c , and larger at p > r c (Gilfanov et al. 1987b). Note 



that this effect takes place already at To < 1, when only 
single resonant scattering is important, and becomes more 
pronounced when To > 1. 

When the optical depth becomes larger than unity, 
multiple resonant scatterings come into play. In the sta- 
tionary scenario, the surface brightness profile of scattered 
AGN emission in a line with To > 1 will be flatter in the 
cluster core compared to a 1/p profile in the t «1 case. 
Also, the total flux of AGN photons scattered toward us 
from the core region will be smaller than calculated in the 
single-scattering approximation by a factor ~ To. However, 
due to the rather similar effects taking place for thermal 
gas emission (see above), the net correction to the ratio 
EW/EWq is expected to be fairly small for t > 1. 

Multiple resonant scatterings also delay the escape of 
AGN photons from the cluster core. As a result, the time- 
dependent scattered surface brightness profiles arising in 
the outburst scenario will be affected. 



4. Numerical simulations 

As a next step, we have performed numerical simula- 
tions for two actual clusters. We used a combination of 
1) a Monte-Carlo code to play the diffusion of photons 
(both those emitted by the AGN and by the gas) through 
multiple reson ant scatterings and 2) the MEKAL code 
(Kaastra 1992 ) as im pleme nted in the software package 
XSPEC vlO (Arnaud 1996 ) to calculate plasma emissiv- 
ity in lines. The line energies and oscillator strengths were 
taken from the lis t of s trong resonance lines of Verner, 
Verner & Ferland ( 1996 ) and the s olar a bundances of el- 
ements from Anders & Grevesse ( |1989| ). The ionization 
fractions were calculated using collisional ionization rate 
fits from Voronov ( 1997), rad iative recombination rates 
from Verner & Ferland ( 1996 ) and diclcctro nic rc combi- 
nation rates from Aldrovandi & Pequignot (1973), Shull 
& van Steenberg fll982|) and Arnaud & Rothenflug (|1985[). 



10 



Sazonov et al.: The past X-ray luminosity of AGN 



Table 1. The brightest X-ray and extreme UV resonance 
lines of the M87/Virgo intracluster gas. 



Ion 


Energy 
(keV) 


Equivalent width 
(eV) 


Optical depth 


Fe XXIV 


0.049 


1 


3.0 


Fe XXIV 


0.065 


2 


6.0 


Fe XXIII 


0.093 


3 


8.5 


VIII 


0.654 


20 


0.6 


Ne X 


1.021 


30 


0.9 


Fe XXIII 


1.129 


10 


3.2 


Fe XXIV 


1.166 


30 


2.7 


Si XIII 


1.865 


10 


1.3 


Si XIV 


2.005 


60 


1.3 


S XV 


2.461 


20 


1.3 


Ar XVII 


3.140 


20 


0.4 


Fe XXV 


6.700 


740 


1.2 


Fe XXV 


7.881 


110 


0.2 



4.1. M87/ 'Virgo 

One of the most promising targets in the sky for observa- 
tions of the discussed effect is M87. This is a giant ellipti- 
cal galaxy located near the center of the nearest (17 Mpc 
- 1 arcmin corresponds to 5 kpc) rich cluster of galaxies 
Virgo. M87 hosts a central black hole of mass 3 x 10 9 M Q 

1997D, which is 



(Macchetto et al. 1997; Marconi et al. 
the largest reliably measured black hole mass. The cor- 
responding Eddington luminosity ^Edd = 4 x 10 47 erg/s. 

M87 demonstrates AGN activity. In particular, obser- 
vations at different energy bands reveal an unresolved nu- 
cleus and a one-sided jet. The combined bolometric lu- 
minosity of the nucleus and jet is ~ 10 42 erg s^ 1 (e.g. 



Biretta, Stern & Harris |1991| ), or ~ 1(T 6 of the Edding- 
ton luminosity. This is an amazingly small value given the 
fact that the nucleus is located in the center of a dense 
cooling flow. The average power input in the form of rel- 
ativistic plasma from the jet into a large-scale radio halo 



around M87 is estimated (Owen, Eilek & Kassim 200C) to 
have been ~ 10 44 erg s _1 over the past ~ 10 s years, which 
is still a small fraction of the Eddington luminosity. 

Both the nucleus and the jet (whose brightest knot is at 
a distance of ~ 1 kpc from the nucleus) have recently been 
observed in X-rays with the XMM-Newton (Bohringer et 
al. |200ip and Chandra (Wilson & Yang |2002| ) satellites. 
The X-ray spectrum is power-law with values for the pho- 
ton index ranging between 2 and 2.9 for the nucleus and 
different knots in the jet. The inferred X-ray (1-10 keV) 
luminosities of the nucleus and the jet are ~ 3 x 10 40 and 
8 x 10 40 erg, respectively (Bohringer et al. [2001 ). 

In our simulations, we modeled the distribution of in- 
tracluster gas around M87 based on the results of recent 
analyses of the XMM observations by Bohringer et al. 
(|2001|) , Finoguenov et al. Q2002[ ) and Matsushita et al. 
( 2002 ). We assumed the distribution of gas to be spher- 



eration to the central region with a radius of 250 kpc. 
The electron density radial profile was represented by a 
sum of two beta-models ([)]) with the following param- 
eters: n — 0.13 cm~ 3 , r c = 1.7 kpc, j3 — 0.42 and 
n = 0.011 cm" 3 , r c = 22 kpc, (3 = 0.47. We adopted 
a constant temperature kT = 1 keV within the central 
1 kpc, the dependence kT = 1.5(r/5 kpc) 0,22 keV be- 
tween 1 and 50 kpc, and a constant kT = 2.5 keV at 
larger radii. The last value approximately correspo nds to 
temperatures measured with ASCA (Shibata et al. 2001 ) 
in different parts of the Virgo cluster within ~ 300 kpc 
from M87. The element abundances were taken to be con- 
stant at r < 10 kpc, with A(0) = 0.4, A(Si) = 1.0, 
A(S) = 1.0, A(Ar) = 1.0 and A(Fe) = 0.8 (in units of so- 
lar abundances), then gradually declining to A(0) — 0.4, 
A(Si) = 0.6, A(S) = 0.6, A(Ar) = 0.6 and A(Fe) = 0.4 at 
r = 40 kpc, and constant thereafter. 

Table [l] lists resonance X-ray and extreme UV lines 
that we estimate to be the strongest emission lines of the 
M87 gas. All of these lines have substantial optical depths 
from the center to the boundary of the gas cloud, assuming 
negligible turbulence. 

Using equation (|20|), we can estimate for the station- 
ary scenario the contribution of scattered AGN radiation 
to the surface brightness of M87/ Virgo in the X-ray con- 
tinuum: 



Dscat 
Be, cont 



Dscat 
-°.E,cont 

Be, cont 



4 x 10 2 



L 



X 



^Edd 



£-°- b exp(£yi.5keV) 



P 



1.7 kpc 



1.7 kpc < p< 22 kpc; (30) 



40^-- £T ' 6 cxp(£;/2 keV) 



-^Edd 
P 



22 kpc 



-0.6 



, p > 22 kpc. 



(31) 



Here we have assumed constant temperatures of kT = 
1.5 keV and kT = 2 keV for equations @ and ©, 
respectively, and 7 = 2. 

The corresponding dependences for the equivalent 
width of a resonance line with to « 1 can be found 
from equations <M) and @ as AEW/EW = (R - 



ically symmetric around M87 and restricted our consid- 



1 )- B l : ,cont(- E, ifc)/- B s,cont(£'ife) (assuming constant element 
abundances for a given p). 

Now we proceed to discussing the results of our compu- 
tations, which do not use the assumption of small optical 
depth. 



4.1.1. Switch-off scenario 

We first consider a scenario in which M87 was bright in 
the past for a long time (at least ~ a million years) until it 
suddenly turned off (actually switched to its present low- 
luminosity state) at t a g — 0. We have plotted in Figs. ||a 
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Fig. 5. a) Ratio of the surface brightness of scattered 
AGN radiation to the intrinsic brightness of the hot gas 
in the iron 1.13 keV resonance line, in the switch-off 
scenario for M87, at different times after AGN switch- 
off as a function of projected radius. The AGN emitted 
Lx = O.OlLEdd in the past. The different solid lines cor- 
respond, from top to bottom, to elapsed times of 5, 50, 100 
and 250 thousand years. The dashed lines are the corre- 
sponding profiles for the spectral continuum at the energy 
of the resonance line, b) Same as Fig. [l]a, but for the iron 
Fe 7.88 keV line, c) Same as (a), but for the iron 6.70 keV 
line and a 10 5 -year outburst (Lx = O.OlLEdd). d) Same 
as (c), but for a 2 x 10 4 -year outburst. 

and |5|b for two resonance lines the computed contribution 
of scattered radiation from the AGN to the surface bright- 
ness as a function of projected radius as would be mea- 
sured by observers living in various epochs t > after the 
switch-off. Also shown are the corresponding profiles for 
the continuum emission near the lines. The relative con- 
tribution of scattered AGN emission is a factor of R ~ 6 
(with a scatter of ~ 3 between the lines) larger for the 
lines than for the continuum, in agreement with the result 
of §0. We adopted Lx = O.OlLEdd for our examples. Since 
the discussed effect is proportional to Lx , the results can 
be easily recomputed for any AGN luminosity. 

Fig. ^ compares the surface brightness profiles of scat- 
tered AGN radiation and thermal emission in the iron 
6.7 keV line, which has the largest equivalent width in 
our sample. It is essential that photon statistics is not ex- 
pected to impose strong restrictions on the prospects of 
detecting scattered AGN radiation from the Virgo cluster 



Fig. 6. M87/ Virgo surface brightness of scattered AGN 
radiation (solid lines) in the iron 6.70 keV resonance line at 
different times (as in Fig. ||) after the AGN switch-off. The 
dashed line shows the intrinsic surface brightness profile 
of the hot gas in the 6.70 keV line, which also is distorted 
by resonant scattering. The AGN emitted Lx — O.OlLEdd 
in the past. 



out to ~ 1 Mpc from M87 with the next generation of 
X-ray spectrometers. 

The finite depth effects that were discussed in §3.3 



play in the M87 case only a minor role for the X-ray lines, 
because for them r < 3 (see Table |l|). Fig. ^ compares 
two cases, M = (no turbulence) and M — 0.5 (signifi- 
cant turbulence), for the iron 1.13 keV line. The increased 
turbulence leads to a decrease from ro = 3.2 to To = 0.7. 
We see that the relative contribution of scattered AGN 
emission to the surface brightness in the central ~ 10 kpc 
is several times larger in the former (optically-thick) case, 
except for the early moments. This increase is mainly due 
to the delayed escape of photons. 

We can now address the principal question of our 
study: what constraints can X-ray observations place on 
the past luminosity of the M87 central source? Suppose 
that it is possible in an experiment to map the distri- 
bution of the equivalent width of a particular line over 
M87/Virgo with 6 — 10% accuracy, which we consider a 
realistic value (see the discussion in §||). Should then the 
relative contribution of scattered AGN emission exceed 6 
at least at some projected radius (i.e. at least at the maxi- 
mum of EW I EWq in Fig. [sj) , it will be possible to estimate 
the past AGN luminosity. The resulting dependence of the 
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Fig. 7. Ratio of the surface brightness of scattered AGN 
radiation to the intrinsic brightness of the hot gas in the 
iron 1.13 keV resonance line, in the switch-off scenario for 
M87, as a function of projected radius. The sampled times 
are as in Fig. |[ The solid lines correspond to the case of 
negligible turbulence, and the dotted lines to M = 0.5. 
Note the linear scaling chosen for the projected radius, as 
opposed to the logarithmic scaling in Fig. ||. 



minimum detectable past luminosity of M87 on time after 
the switch-off is plotted in Fig. [8|a. 

We see that in order to have a 10% contribution of 
scattered radiation from the AGN whose luminosity in 
the past was L x = 10~ 4 LEdd = 4 x 10 43 erg s _1 (i.e. 
several hundred times higher than the combined X-ray 
luminosity of the nucleus and jet now), we should be not 
later than ~ 3x 10 4 years after the switch-off and look with 
X-ray telescopes into the M87 core region (p < 10 kpc ~ 
2 arcmin in the plane of the sky), as is clear from Fig. If 
the luminosity in the past was higher, Lx = 10~ 3 LEdd = 
4x 10 44 erg s _1 , we must live not later than ~ 2x 10 5 years 
after the switch-off and look further out (at p ~ 50 — 
200 kpc ~ 20 arcmin from the AGN). The nearly flat 
(oc t~ - 6 ) shape of the curves in Fig. |]a, which have been 
truncated at t = 3 x 10 5 years to mimimize the effects of 
the finite size of our model gas cloud, suggests that if the 
X-ray luminosity of M87 was ~ a few 10 _3 £Edd ~ 10 45 erg 
s~ 4 in the past, the maximum allowed elapsed time would 
only be restricted by the extent of the gas in the Virgo 
cluster (~ 1 Mpc), the effective area of the detector and 
the field of view of the telescope. 



Fig. 8. a) Minimum detectable (for S = 10% - see 
text) past X-ray luminosity as a function of elapsed 
time in the switch-off scenario for M87. The different 
curves correspond to mapping M87/Virgo in different 
lines: Fe 1.13 keV line (solid), Fe 1.17 keV line (dotted), 
Si 1.86 keV line (short-dashed), S 2.46 keV line (long- 
dashed), Ar 3.14 keV line (dot-short-dashed), Fe 6.70 keV 
line (dot-long-dashcd) and Fe 7.88 keV line (short-dash- 
long-dashed), b) Same as (a), but for a 10 5 -year flare, c) 
Same as (b), but for a 2 x 10 4 -year flare. 



We may therefore conclude that it should be possi- 
ble with future X-ray spectrometers to measure or obtain 
fairly tight limits on the past X-ray luminosity of the cen- 
tral source in M87. For comparison, we can use the avail- 
able data of XMM-Newton observations of the M87/Virgo 
core to derive some weaker limits already now. Bohringer 
et al. (2001) have analyzed X-ray spectra taken from a 
sequence of rings centered on the M87 nucleus. The most 
distant ring has an outer radius r ou t ~ 60 kpc, and the 
innermost circle has r- m ss 1.25 kpc. Except for the spec- 
trum of the central circle to which the direct emission from 
the AGN contributes, the individual spectra correspond to 
optically thin thermal emission from the gas. The quality 
of the spectra is apparently good up to ~ 7 keV. We can 
therefore suggest that the surface brightness of scattered 
AGN radiation (having a power-law spectrum) does not 
exceed the actually measured surface brightness in the X- 
ray continuum at 7 keV nowhere between p = 1.25 kpc 
and 60 kpc. From this condition we can find an upper 
limit on the past X-ray luminosity of M87 as a function of 
time after the switch-off, similarly as we derived the limits 
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Table 2. The brightest resonance X-ray and extreme UV 
lines of the Cyg A intracluster gas. 




50 100 150 200 250 

Elapsed time (10 3 years) 



Fig. 9. Upper limits (for 5 = 100%) on the past X-ray 
luminosity of the M87 AGN as a function of elapsed 
time, derived from XMM-Newton observations of the 
M87/Virgo central region in the X-ray spectral contin- 
uum below 7 keV, for three cases: switch-off, a 10 5 -year 
flare and a 2 x 10 4 -year flare (solid lines from bottom to 
top). For comparison shown (dashed lines) are the lim- 
its obtainable if it were possible to map M87/ Virgo at 
E = 20 keV. 



presented in Fig. ga based on the resonance line compu- 
tations. However, in the previous case we did not impose 
any limitations on the projected radius and we adopted 
a value 5 — 10%, rather than 5 = 100% to be used now. 
The resulting dependence is shown in Fig. |^. 

We see that the current limits based on the spectral 
continuum data are rather weak. For example, the aver- 
age X-ray luminosity of the M87 AGN over a long period 
of ;> 10 5 years ending 5 x 10 4 years ago was not higher 
than 5 x 10~ 3 LEdd = 2 x 10 45 erg s — 1 . By comparing 
Fig. | with Fig. | a, we see that the current limits could 
be improved by at least an order of magnitude by means of 
fine X-ray spectroscopy. Alternatively, much tighter limits 
could be obtained by measuring the surface brightness of 
M87/ Virgo in the hard X-ray continuum. To illustrate this 
possibility, we have repeated our analysis for E = 20 keV 
- see the result in Fig. ||. Going from E = 7 keV to 20 keV 
(again assuming a S — 100% scattered contribution rela- 
tive to the bremsstralung emission) leads to a factor of 
~ 200 (!) stronger limits. 
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4.1.2. Flare scenario 

Suppose now that some time ago M87 experienced an out- 
burst that was short compared to the characteristic light 
crossing time of the gas cloud, i.e. A <C 10 6 years. Some 



analytic estimates for this case were obtained in § 3.2.3 
Figs. |c and |d show the computed time evolution of 
the scattered/thermal brightness ratio profile for the Fe 
XXV Ka line for outburst durations A = 10 5 years and 
A = 2 x 10 years, respectively, and an AGN luminos- 
ity Lx = O.OliEdd- Note that for t ^> A the scattered 
surface brightness profile depends on the product Ljc A, 
rather than on Lx and A separately. 

Figs. ||b and ||c show the minimum detectable (for 
10% scattered contribution) AGN outburst luminosity as 
a function of elapsed time for the two considered scenar- 
ios. These dependences should be compared with Fig. |^a 
representing the switch-off scenario. We see that until a 
certain moment t cr it ~ A the minimum detectable lumi- 
nosity is only a little smaller in the flare case than in 
the switch-off case. This is expected, because the char- 
acteristic line-of-sight depth of the illuminated volume at 
t = £ C rit is about the same, bp ~ ct cr j t , in both cases. Only 
at t 3> A, does the minimum detectable luminosity be- 
come much larger (by a factor of A/i) for short outbursts. 
This directly follows from comparison of equations (20) 
and (||). 

Fig. ^ shows upper limits for the same outburst sce- 
narios but obtained using the XMM-Newton data on the 
X-ray continuum emission from M87/Virgo. 

4-2. Cygnus A 

Our second example is Cyg A. This is a well-known nearby 
(z — 0.0562, i.e. about 20 times more distant than M87 - 
1 arcmin approximately corresponds to 100 kpc) powerful 
radio galaxy. Cyg A has recently been observed with the 
Chandra satellite, and detailed information was obtained 
on the morphology of X-ray emission from the galaxy 
(Young et al. |2002| ) as well on the intracluster gas sur- 
rounding it (Smith et al. 2002). 

In the Chandra an d qua si-simultaneous RXTE obser- 
vations (Young et al. 2002 ), hard (up to 100 keV), spa- 
tially unresolved X-ray emission was detected from the 
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position of the radio and infrared nucleus of Cyg A. The 
energy spectrum of this radiation is power law with pho- 
ton index 7 = 1.5, heavily absorbed below a few keV. The 
inferred unabsorbed 1-10 keV luminosity of the nucleus 
Lx ~ 5 x 10 44 erg s _1 . This is four orders of magnitude 
more than the X-ray luminosity of the M87 nucleus. 

The Chandra observations clearly reveal an intraclus- 
ter medium. The gas has complex structure within the 
central ~ 100 kpc, which apparently is the result of an in- 
teraction with the relativistic material produced by AGN 
activity. In contrast, the morphology of the gas is simple 
at larger radii (;> 100 kpc), namely the gas appears to be 
spherical within at least 700 kpc of the nucleus, and nearly 
isothermal with kT ^7 — 8 keV. The total luminosity of 
the intracluster gas is ~ 10 45 erg s _1 , which is compara- 
ble to the X-ray luminosity of the central AGN and is a 
typical value for rich clusters of galaxies. 

Given the above observational facts, we can estimate 
the contribution of scattered emission from the AGN to 
the X-ray surface brightness outside the central 100 kpc of 
the Cyg A cluster. Due to the highly irregular distribution 
of gas in the innermost region of Cyg A, our simulations 
will not be aimed at this zone. Based on these considera- 



tions and the results of Smith et al. ( |2002D , we model the 
radial distribution of gas density by a beta-model with 
no = 0.05 cm -3 , r c = 30 kpc and f3 = 0.5. The tem- 
perature is assumed to be constant, kT — 5 keV, within 
the central 50 kpc and also at r > 150 kpc, kT = 7.5 keV. 
The abundance of iron is taken to be 0.35 solar. Our model 
cloud of gas has an outer boundary at r = 1 Mpc. 

Table || lists the strongest emission lines of the Cyg A 
cluster, as implied by our model. The Chandra and RXTE 
data do reveal two strong lines, one near 6.7 keV (Fe Ka) 
and another near 7.9 keV (Fe K/3 pl us po ssibly Ni Ka). 
It appears (see Fig. 7 in Smith et al. |2002| ) that the mea- 
sured (with moderate energy resolution) spectra also do 
not contradict the presence of a strong line at the position 
of Fe Lya (6.97 keV). 

We can estimate from equation (^o|) the stationary 
contribution of scattered AGN radiation to the surface 
brightness of the Cyg A cluster in the X-ray continuum: 



Dscat 
-°,E,cont 



D 



2 X 10" 



L 



X 



S.cont 



5 x 10 44 ergs" 
-0.5 



30 kpc 



J B" ul exp(£;/7.5keV) 



(32) 



where we have normalized the X-ray luminosity of the 
Cyg A nucleus to its present-day value and used 7 = 1.5. 

Fig. [lO] shows a computed Cyg A radial profile of 
the scattered/thermal surface brightness ratio for the iron 
6.7 keV line, assuming that the AGN X-ray luminosity has 
remained for a few million years the same as it is now (the 
stationary scenario) . We see that scattered AGN emission 
contributes of the order of 1% at p ~ 200 kpc. It then fol- 
lows that the Cyg A luminosity should have been an order 
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Fig. 10. Ratio of the surface brightness of scattered AGN 
radiation to the intrinsic brightness of the hot gas in the 
iron 6.70 keV resonance line, in the stationary scenario for 
Cyg A, as a function of projected radius (the solid line). 
The AGN X-ray luminosity is assumed to be equal to its 
present-day value, Lx — 5 x 10 44 erg s _1 . The dashed line 
shows the corresponding profile for the spectral continuum 
at the resonance energy. The vertical dotted line indicates 
that the results to the left of it are not reliable. 



of magnitude higher in the past (i.e. Lx ~ 5 x 10 45 erg 
s _1 ) for scattered AGN radiation to contribute ~ 10% to 
the strong X-ray lines observed from the intracluster gas 
(see Fig. [TTj> . 

We can also use the available Chandra data on the X- 
ray continuum emission from the Cyg A cluster to find 
upper limits on the past X-ray luminosity of its cen- 
tral source. We adopt the following parameters for this 
analysis: E = 8 keV - the maximum photon energy at 
which thermal bremsstrahlung emission is still detected, 
and p m in = 10 2 kpc and p max = 6 x 10 2 kpc - the bound- 
aries of the region that is well explored with Chandra. We 
plot the resulting upper limits in Fig. [Tl|b. We see that 
the current constraints are very weak because of the high 
temperature (~ 7 keV) of the intracluster gas. For com- 
parison (see Fig. p"T|b), arcminute-resolution observations 
of the Cyg A cluster at E ~ 40 keV could provide lim- 
its on its past luminosity similar to those obtainable by 
mapping the cluster in resonance lines (see Fig. |Tl|a). 
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Fig. 11. a) Minimum detectable (for 5 = 10% - see 
text) past X-ray luminosity of Cyg A as a function of 
elapsed time. The surface brightness is measured in the 
iron 6.70 keV line. The solid line corresponds to the 
switch-off case (A — > oo). The dotted and dashed lines 
correspond to outbursts of duration A = 5 x 10 5 years 
and A = 2 x 10 5 years, respectively. The luminosity is nor- 
malized to the present-day X-ray luminosity of the Cyg A 
nucleus. The vertical dotted line indicates that the results 
to the left of it are not reliable, b) Same as (a), but de- 
rived from Chandra observations of the Cyg A cluster in 
the X-ray spectral continuum below 8 keV (for 5 = 100%). 
For comparison shown (the lower set of lines) are the lim- 
its that could be obtained if it were possible to map the 
Cyg A cluster at E = 40 keV. 



5. Discussion and conclusions 

The main results of this paper are as follows. 

— We have demonstrated that there are two viable ob- 
servational strategies for constraining the past X-ray 
luminosity of galactic nuclei with the next generation 
of X-ray telescopes. One is to search for scattered AGN 
radiation in the spectral continuum at high energies 
E ^> kT(l + z), and the other is to perform imaging 
in bright resonance X-ray lines. The first approach ap- 
pears to be particularly promising for studying distant 
(at z > 0.5) powerful quasars and their environments 

(§ED- 

— We have shown that the relative contribution of scat- 
tered radiation is typically larger by a significant fac- 



tor of 3-10 in a resonance line than in the neighboring 
continuum (§Q). 
— We have estimated the level of constraints that could 
be derived from future observations on the past X-ray 
luminosity of the nearby M87 and Cyg A active galax- 
ies (§||). The upper limits imposed by the available 
XMM-Newton and Chandra X-ray continuum data are 
typically 1-2 orders of magnitude weaker. 

Although we have been mostly discussing clusters of 
galaxies and their central dominant galaxies, the same 
method can also be applied to groups of galaxies and iso- 
lated giant elliptical galaxies, which also are large reser- 
voirs of hot ionized gas capable of scattering AGN radia- 
tion. The main advantage of clusters is the larger extent of 
intracluster gas, which enables studying previous galactic 
activity on longer timescales (up to a few times 10 6 years 
for the richest clusters) than in the case of groups of galax- 
ies and early- type galaxies (up to a few times 10 5 years). 
On the other hand, the latter two classes of object have 
their own advantage that typical gas temperatures for 
them are kT ~ 1 keV and the factor exp(E / kT) becomes 
significant already in the standard X-ray band (< 10 keV). 
In particular, one can use the Ka line of He-like iron at 
6.7 keV. Indeed, the corresponding ions (Fe XXV) are still 
abundant at such low temperatures, but collisional exci- 
tation of the ions leading to emission in the line is already 
inefficient. 

Although we have only considered the case where the 
AGN is an isotropic source, in reality angular anisotropy 
at some level is expected. Therefore, the contribution of 
scattered emission from the AGN to the surface brightness 
may be larger or smaller in some regions than predicted as- 
suming source isot ropy. W e refer the reader to the p apers 
by Gilfanov et al. fll987aj ) and Wise & Sarazin fll992| ) for a 
detailed discussion of scattered surface brightness profiles 
that may arise in the case of beamed source emission. 

In the discussion of our simulations for the M87 and 
Cyg A galaxies and their associated clusters of galaxies we 
used a fiducial value of 10% as the minimum detectable 
contribution of scattered AGN emission to a resonance 
line. We realize that significant efforts will be required to 
achieve this or better detection level experimentally. One 
difficulty is that the measured surface brightness is the in- 
tegral of emission (thermal plus scattered) along a given 
line of sight, and there may be present gas with differ- 
ent temperatures and element abundances along this line 
of sight. Both types of variations will have an effect on 
the equivalent line width. In principle it should be pos- 
sible to determine and subtract the underlying tempera- 
ture and abundance radial trends from a careful analysis 
of the spectral-imaging data. The remaining uncertain- 
ties are expected to be largest for cluster central regions, 
where both the gas temperature and element abundances 
vary significantly, for example in the M87 case. On the 
other hand, this problem is not expected to be severe 
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outside cluster cores. For example, in the Cyg A cluster 
the gas appears to be nearly isotropic outside the central 



~ 200 kpc (see Smith et al. 2002), whereas the equiva- 
lent width of the 6.7 keV iron emission line changes by 
only ~ 10% when the gas temperature varies by 0.5 keV 
around kT ~ 7.5 keV. Therefore, achieving the 10% de- 
tection level appears realistic. 

It is clear that the potential of the method can be fully 
realized only when it is possible to resolve the interest- 
ing resonance lines from neighbouring lines. Indeed, lines 
other than resonance lines can contribute significantly to 
the intrinsic emission of the gas, but little to the scattered 
AGN emission. Moreover, it is desirable that all important 
lines in blends be resolved, such as the permitted, inter- 
combination and forbidden lines in the complex around 
the 6.7 keV (Fe XXV Ka) line. Even a small contribution 
of scattered emission from the AGN could then be made 
manifest by comparing the surface brightness profiles mea- 
sured in these lines. We are therefore looking forward to 
future high-energy astrophysics missions such as Astro-E2, 
Constellation-X and XEUS that will provide the required 
energy resolution (~ a few eV). 

Finally, we would like to note that the constraints on 
the past AGN X-ray luminosity could be further improved 
by means of X-ray polarimetry, especially if it became pos- 
sible to measure polarization in X-ray resonance lines. The 
radiation of a central source gets strongly polarized upon 
scattering in a beta-cluster (with a typical resulting de- 
gree of polarization P ~ 60%, Sunyaev 1982), whereas 
the intrinsic emission of the intracluster gas in X-ray res- 
onance lines is also polarized but to a much lesser degree 
(P <; 10%, Sazonov et al. 2002). Thus, even a few per 
cent contribution of scattered AGN radiation to the sur- 
face brightness in a resonance line would be manifest in 
polarimetric observations. 
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